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Abstract 



. . -This paper analyzes coherent backscattering of intense laser light by two randomly placed distant atoms. Starting from the general 
two-atom master equation, we analytically derive the elastic and inelastic background and interference components of the double 
scattering spectrum. By expressing the final results in terms of single-atom observables, the two-atom problem is shown to be 
equivalent to a description in terms of single atoms under bichromatic driving. 
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^T) 1. Introduction 

I— !■ Coherent backscattering (CBS) of light is a phenomenon 
emerging in a dilute disordered medium, whereupon multiply 
II scattered counter-propagating waves interfere constructively in 
a narrow angular range around the backward direction, lead- 
ing to enhanced intensity of the scattered signal fT]. CBS of 
^ light has attracted ongoing interest due to its close connection 
O^to various aspects of interference-induced effects in presence of 
disorder, such as weak localization [5], Anderson localization 
^— I or random lasers 

' Since 1999, CBS of light has been observed in cold atoms 
Jsl]. Quantum scatterers are endowed with internal structure 
, that impacts the matter-light interaction already at the level of 
\^ a single atom and, consequently, influences multiple scatter- 
ing. On the one hand, in the elastic scattering regime of low 
laser intensity, the degeneracy of atomic dipole transitions can 
^trongly reduce the CBS signal, and needs to be taken into ac- 
count into the theory in order to meet the experiment |3, 01 • 
^ On the other hand, for degenerate (e.g., in rubidium atoms) 
• ^ and non-degenerate (e.g. in strontium atoms ^) transitions 
alike, increasing the intensity of the incident laser field leads 
^ to atomic saturation accompanied by nonlinear (multi-photon) 
" " inelastic scattering processes ^ reducing the interference con- 
trast in CBS (id]. 

While a thorough theory of CBS of light from de gen erate 



m 
o 
o 



atoms in the regime of elastic scattering already exists 111 Ull2ll . 



a theory of CBS from saturated atoms is still in a rudimen- 
tary state, several years after the "saturation" experiments with 
strontium ifioll and rubidium [13] atoms. The main obstacle for 
developing such a theory is that CBS becomes a multi-wave 
interference effect in the nonlinear scattering regime [ 14, ,15,1 . 
with more than two interfering amplitudes. A nonperturbative 
(with respect to the strength of the nonlinearity) theory of non- 
linear CBS that fully accounts for the multi-wave interference 
character, has so far been developed only for nonlinear classi- 
cal scatterers lfl4l[l5ll . not including quantum effects resulting 



from multi-photon scattering processes in the regime of atomic 
saturation. 

Among the present theoretical approaches to CBS from sat- 
urated atoms, one lfl6ll is based on diagrammatic scattering the- 
ory and, hence, viable to treating the multiple scattering pro- 
cesses of arbitrary order It is, however, bound to a weakly in- 
elastic regime, and, so far, has been worked out for two incident 
laser photons (l^]. Other approaches are based on quantum op- 
tical master [17^ 20.1 or Langevin lilSil equations, and treat in 
detail the response of individual atomic scatterers to the inci- 
dent field of arbitrary intensity, though are restricted to a few 
atoms. 

Recently, we have proposed a new approach to multiple scat- 
tering of light from distant atoms in the saturation regime 12 ill . 
This method, to which we shall refer to as 'the pump-probe' 
approach, unifies the potential of diagrammatic scattering the- 
ory with that of quantum optical methods, and, as we believe, 
will enable us to develop a multiple scattering theory of intense 
laser light from dilute atomic gases, where the distances be- 
tween atoms are much larger than the laser wavelength. 

Within the pump-probe approach, which has so far been elab- 
orated in detail for two scalar atoms ll2llE2 [l. the configuration- 
averaged double scattering signal is extracted from single-atom 
observables (such as the spectral correlation function of the 
atomic Bloch vector), where each atom is subject to a bichro- 
matic driving field, consisting of the laser field (pump) and the 
weak field scattered from the other atom (probe) whose fre- 
quency may differ from the laser frequency. The main assump- 
tion of the pump-probe approach is that both these fields can 
be modelled as classical fields. Intuitively, this is expected to 
be possible in the case of large distance between both atoms, 
where only a single photon is scattered from one atom to the 
other, and thus correlations between different photons emitted 
by the same atom play no role. 



As the numerical data shows 112 iL 12211 . the double scatter- 



ing spectra calculated within this new method are in excellent 
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agreement with the accurate results following from the two- 
atom master equation expanded up to lowest nonvanishing or- 
der in the inverse distance between both atoms. This shows 
that the ansatz of a classical probe field is indeed justified for 
large distances. Furthermore, for the elastic component of the 
backscattered field, we have analytically proven the strict equiv- 
alence of these two apparently very different approaches 12211 . 
However, the full analytical proof including also the inelastic 
spectra has been missing. 

In the present contribution, we fill in this gap. Starting from 
the master equation for two scalar atoms, we derive the double 
scattering CBS elastic and inelastic signals, and express them 
in terms of purely single-atom quantities. Next, we establish 
the equivalence with the result of the pump-probe approach, 
based on single-atom equations under the bichromatic driving 
|2lil22|l. 

The structure of the paper is as follows. In the next section, 
we introduce the master equation approach and define all quan- 
tities of physical interest. In Sees. [3]and|4]we deduce, respec- 
tively, the elastic and inelastic components of the backscattered 
light from two atoms. In Sec.|5]we interpret the obtained results 
by considering a physical setting including single atoms sub- 
jected to bichromatic driving. Our work is concluded in Sec. |6l 

2. Master equation approach 

2.1. Formalism 

Let us consider a toy model of CBS consisting of two scalar 
two-level atoms embedded in a common electromagnetic bath 
and driven by a quasi-resonant laser wave with a wave vector 
ki- In the following, we will neglect the atomic center-of-mass 
motion and focus on the dynamics of the internal degrees of 
freedom. The standard method to find the atomic dynamics is a 
master equation approach ll23il24tl . whereupon the bath degrees 
of freedom are traced out leading to the equation of motion for 
the quantum mechanical expectation values of atomic variables. 

It should be noted that a more general model taking into ac- 
count the vectorial nature of the scattered field as well as of 
the atomic transitions has already been studied in some detail 
l25ll . The reason why we are addressing here the simpler 
scalar case is that we want to analytically derive expressions for 
the power spectrum of light coherently backscattered from two 
atoms, and compare them with the results of 1I21I1 obtained for 
scalar atoms. However, we stress that the tools developed in 
this work can be generalized for arbitrary atomic transitions. 

In the Heisenberg picture and in the frame rotating at the laser 
frequency, one obtains the following master equation for the 
expectation value of an arbitrary atomic operator Q ll23ll : 



Here, crj = |1);(2|^' and crt = |2>y<l|^', with and \2)j be- 
ing respectively the ground and excited states of atom j, de- 
note the atomic lowering and raising operators. Furthermore, 
D.j - Q.e'^'- ''' is the Rabi frequency dependent on the atomic 
position rj, S - - uq \s the laser-atom detuning, and y is 
half the Einstein's A coefficient. The lower line of Eq. ([TJ de- 
scribes the retarded dipole-dipole interaction dependent on the 
dimensionless parameter jc = (jJ{)\r\ - ral/c. Coherent backscat- 
tering can be observed in the dilute regime x » 1, where the 
atoms are located in the far field of each other and exchange a 
single photon. Within our toy model, this corresponds to a per- 
turbative treatment of the dipole-dipole interaction up to second 
order in the interatomic coupling T{x). The explicit form of the 
latter in the limit jc » 1 can be presented as 



3 e-" 

T{x) ^ly- . 
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(2) 



By inserting into Eq. ([T]) operators from the complete two- 
atom basis set of operators: Q ^ {q\® qi], with 



h 



(3) 



one obtains a closed linear system of 16 coupled equations for 
the atomic averages. Since {![ (8/2) = 1 for every atomic state 
(due to normalization), this can be reduced to a 15 dimensional 
system for the two-atom Bloch vector 



with 



<0 = «<?2),(<?l>,<<?l®<?2>y 



(4) 



(5) 



This system of equations has the following matrix representa- 
tion: 

<^) = (A + y)<^) + L+,o, (6) 

where the matrices A and V describe the evolution of indepen- 
dent and dipole-dipole interacting atoms, respectively: A (as 
well as the free vector L+_o) is generated by the two upper lines, 
and V by the lower line of equation ([T]|. The explicit forms of 
the matrices A and V as well as vector L+ will be given in 
SeclH 



2.2. Intensity of coherently backscattered light 

The solution of Eq. (|6]l not only yields the atomic evolution, 
but also the field radiated by the atoms. For example, the av- 
erage stationary scattered light intensity in the direction of the 
unit vector n. 



-y{cr}crjQ + Qcr-crj-2cr^Qcrj)} 

2 



(1) 



(m) = lim{E(-\t)E^^\t)h. 



(7) 



after omitting unimportant prefactors and neglecting retarda- 
tion, is given by 

m) = Y.^cr^'^jhs + {^}T,)sse''"'<''-'^\ (8) 
J JMj 



where the subscript means that the expectation value is evalu- 
ated in the steady state. Bearing in mind that, in the following, 
all quantum mechanical expectation values are to be evaluated 
in the steady state, we will henceforth lighten notation by drop- 
ping this subscript. 

We note that the light intensity, Eq. (O, radiated by two 
distant atoms has already been considered previously, for ex- 
ample in 1 26] (where the two atoms are coupled by an addi- 
tional optical element). However, in the following, our focus 
will not be on the total intensity, as given by Eq. ([8]l, but, with 
view at the multiple scattering scenario from a dilute cloud of 
atoms, we are rather interested in the double scattering coher- 
ent backscattering intensity, which is extracted from Eq. ([8]l by 
expansion in second order in the atom-atom interaction and sub- 
sequent average over the positions of atoms. Zeroth-, first-, and 
second-order terms of the perturbative expansion of Eq. Q in 
the coupling Y will carry upper indices '(0)', '(1)\ and '(2)', 
respectively. That is, all quantities that will be studied here are 
quadratic in the dipole-dipole coupling, and will be supplied 
with superscript '(2)'- For example, the double scattering in- 
tensity is obtained from 

</(„^)(2) . Yia)a-,r + 2<->A->®^"'"''"''''- (9) 

i iMj 

Eq. ^ does not yet correspond to the CBS intensity. The latter 
is obtained after the configuration averaging procedure, denoted 

as (. . .)conf., 

«/(^)<'')conf. = Yj^{(r]cT-^f\,nf. (10) 

j 

JMj 

and includes two terms: the background, or ladder intensity, 
given by the first term in the right hand side of Eq. ( fTOl i. and 
the interference, or crossed intensity, given by the second part 
of Eq. ([Tol l. The background and interference CBS intensi- 
ties in Eq. (fTOb are formed respectively by co- and counter- 
propagating double scattering amplitudes, cf. Fig.|2]below. The 
background intensity is radiated uniformly into all directions n, 
whereas the interference part contributes only in the backscat- 
tering direction, with a small angular width AO cc 1/x <s 1. 
Thus the full angle dependence of the coherent backscatter- 
ing intensity is determined by the background intensity and the 
height of the interference peak in exact backscattering direc- 
tion, i.e. n parallel to -ki, which we will assume from now 
on. 

The configuration averaging procedure includes two stages. 
At one stage, averaging over interatomic distances is per- 
formed. As a result, the oscillating terms oc T^(x),(T*(x))^ 
will vanish, while terms oc |r(x)p which vary smoothly with x, 
will be preserved. Another stage of configuration averaging in- 
cludes averaging over random orientations of the radius-vector 
connecting the atoms (to kill oscillations of Qj, Q* and most of 
their combinations). Note, however, that the phase of oscillat- 
ing terms oc fi^Q^ cancels itself with that of the phase factor 
^ikftifj-i'i,) jjjg backwards direction (second term of Eq. Q), 



since k x ki. This is precisely the interference CBS contribu- 
tion surviving the disorder averaging. It will be shown below 
that the averaging over both, distance and orientation, can be 
performed analytically within the master equation approach. 

Once the average CBS intensity is known, one can assess the 
main measure of phase coherence called the enhancement fac- 
tor |6], defined as the interference divided by the background 
intensity, where, as mentioned above, the interference intensity 
is evaluated in exact backward direction, where it assumes its 
maximum value. In the presence of inelastically scattered pho- 
tons, a more refined measure of phase coherence is provided 
by the spectrum of backscattered light |25]. Indeed, the back- 
ground and interference components of the CBS spectrum in- 
dicate the interference character and magnitude as a function 
of the frequency of the backscattered photon; the enhancement 
factor follows from them after integrating over the whole fre- 
quency distribution. 

2.3. Spectrum of coherently backscattered light 

The master equation ([T]l is obviously Markovian, conse- 
quently, due to the quantum regression theorem |9], the multi- 
time atomic correlation functions obey the same equation as 
Eq. (|6]l but with modified initial conditions and free vector. 
Among these multi-time correlation functions, we will consider 
the stationary first-order temporal coherence function of scatter- 
ing atomic dipoles: 

ri(T) = ^«cr;(r,-(r))(2)e'*0'),„„f., (1 1) 

where k = kin, and rjk = rj - rk- 

2.3.1. Elastic spectrum 

The factorized correlation function of atomic dipoles 

r,ei(T) = Y.{i<^T(^k(r)f' + <cr;)<')(^,(r))<" 

+<cr;)(2)<^^-(r))m)e'*''^')^^^^, (12) 

gives the elastic component of the spectrum (oc dico - cjl)) 
through the Laplace transform |0] 

Re r°° 

Sei(co)= fl'Te'"Ti;el(T), (13) 

^ Jo 

since the temporal evolution of the average dipole lowering op- 
erators in the steady state reduces to e^"^'-'^ . Then 

S,x(u>)^lf^6(u>-cjL), (14) 

with 
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The expression inside (. . .)conf. explicitly depends on the co- 
ordinates of atoms j and k. However, after the configuration 
averaging, the result must be entirely symmetric with respect 
to indices interchange j ^ k. Hence, instead of Eq. (fTSl l. we 
can use shortened expressions. Let us write these expressions 
down separately for the background (Lei) and interference (Cei) 
contributions, assuming exactly the backward direction of ob- 
servation (k = -ki): 



Lei - 2{{crir{a-^r> + {(r^.fHcr^f^ 

+ {(rif\cT^f\oni., (16) 

Cei = 2<«^r)^°^K)*'' + <'^r)*"K)<" 

+ <crr)<2'<(r^)(»')e*-'^-)eo„f.. (17) 

Note that in Eq. ( fTSI l. we have fixed j - k - 2, while in Eq. (fTTI l. 
j - 2 and k - 1. 

2.3.2. Inelastic spectrum 

In order to evaluate the inelastic spectrum of the backscat- 
tered light, let us introduce the time-dependent vector 



(18) 



where s'/t) = {cr^QiT)). According to the quantum regression 
theorem the vector As, obeys the equation 



As,- = (A -H y)As,- 
with the initial condition 

As'/O) = <cr;^> - <cr}><^). 



(19) 



(20) 



The inelastic spectrum follows from the fluctuating part of the 
correlation function, Eq. (fTTT i: 



Ari;i„el(T) =ri(T)-ri;el(T) 



= ([As-*f>(T)]4 + [Asf\T)]i 
+ [A^l^\T)he-''-'-- + [A^j'\T)Ue''^-H 



(21) 



COllf. 



via Fourier transformation 



S mel( 



^ Jo 



dre-'^Ari.iM, 



(22) 



where v = co-lji denotes the detuning with respect to the laser 
frequency, and the fact that Eq. ^ is written in the rotating 
frame is taken into account. The indices [. . . ]i and [. . . ]4 in 
Eq. in} refer to the components cr^ and cr^ of the 15 dimen- 
sional vector ^, see Eq. (|4]i. 

Employing the same symmetry argument as we did when 
writing down the expression for the elastic contributions, let 
us write down the corresponding expressions for the inelastic 
background and interference spectra: 



Linel(v)=— I fl'Te'''^<[Asf>(T)],>, 



^ Jo 



conf. 



(23) 



Cinel(v) 



_ 2 Re r 
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fl'Te"'^([As-t,^'^(T)]4e'*^-^'^>eo„f. (24) 



Eqs. ( l23T l and (l24l i can be cast into a more convenient form by 
using Laplace transforms (which exist also for purely imaginary 
Laplace transform variable z - -iv because, as a consequence 
of the fluctuation-dissipation theorem [27] . the fluctuating parts 
of the atomic correlation functions are exponentially decaying 
with time). One obtains 



ii„el(v)= ^<[As7>(-/v)]l), 



conf.? 



Cinel(v) 



2 Re 



{[AS^'\-iv)Ue'^-H 



conf. » 



(25) 
(26) 



where i(z) - J^^ dtexp(-zt)x(t) indicates the Laplace image of 
x(t). 

It follows from Eqs. (O, ([17] (|25]) and (|26]l that finding the 
elastic and inelastic spectra of CBS requires the solution of 
Eqs. (|6]l and (fT9] l. expanded up to second order in the dipole- 
dipole interaction matrix V, with subsequent configuration av- 
eraging of the result. Both of these tasks can be accomplished 
by exploiting the tensor structure of matrices A and V to be dis- 
cussed in the next section. 

2.4. Structure of the evolution matrices 

Let us now study the structure of the matrices A and V intro- 
duced in Eq. (|6]l, referring to the evolution of non-interacting 
and interacting atoms, respectively. 

2.4.1. Free evolution: matrix A 

With the ordering of the 15-dimensional vector defined by 
Eq. (HI, the matrix A receives the following block structure 



A = 



Lx 







(27) 



The block matrix M+ has dimensions 6 x 6; it describes the 
individual evolution of atoms 2 and 1 independently of each 
other: 

M+^M2®Mu (28) 

with the matrix Mj being the standard Bloch matrix for a single 
atom, generated by the upper two lines of Eq. ([T]i: 



-y + i6 -/Q^/2 ^ 
-y-i6 iQ.*/2 



-iO.* 

J 



iQ.i 



-2y 



(29) 



We remind that the optical Bloch equation for a single atom 
reads [28] 



= Mjid'j) + L, L = (0, 0, -2yf. 



(30) 



The 9x9 block Mx describes the evolution of two-atom corre- 
lation functions of uncoupled atoms and, hence, reads 



Mx = Mi ® /2 H- /i ® M2 



(31) 



where Ij is the 3x3 unit matrix. In the following, we will drop 
the indices of the matrices unless the latter diff'er for atoms 1 
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Figure 1: Diagrammatic representation of tlie elementary dipole-dipole in- 
teraction processes between two atoms (black spots correspond to their elec- 
tronic states after the interaction) with the respective terms of the matrix V, see 
Eq. )34) . Solid (dashed) arrows depict positive (negative) frequency photons. 
(a),(d): excitation is lost (gained) by atom 1 (2); (b),(c): same in reversed order. 



and 2. Accordingly, the non-diagonal 9x6 block matrix Lx can 
be written as 

Lx = (L^I I®L), (32) 

with L defined in Eq. ( l30l l. Finally, it is appropriate to specify 
here the vector L+ o introduced in Eq. 



L+,0 = (L^,L^,0,...,0)^. 

9 zeroes 



(33) 



To conclude, we have defined all ingredients of the composite 
system of two non-interacting atoms in terms of single-atom 
Bloch vectors and their tensor products. We will next do the 
same for the interacting atoms. 

2.4.2. Dipole-dipole interaction: matrix V 

We remind that the matrix V is generated by the third line of 
Eq. ([T]i. It naturally splits into 4 components: 



y = yi2 + y2i + Vi*2 + V21, 

'^T(x) ccT'(x) 



(34) 



where the first and second indices in the subscripts coincide, re- 
spectively, with the values of the indices /' and k in Eq. ([TJ. The 
four matrices V12, Vi*2, V21 and V^j describe four elementary 
exchange excitation processes between the atoms which are de- 
picted diagrammatically in Fig. [T] Furthermore, we associate 
terms oc T*{x) (or T{x)) with positive (or negative) frequency 
photons, in analogy with the single-atom case, see Eq. ( l30t . 
where (crj) (or (cj)) are coupled to (cr.) by a coefficient pro- 
portional to the positive (or negative) frequency component of 
the laser field Qj (or D.*). Likewise, in the present two-atom 
case, it is easy to show that the expectation values (crp (or 
(crp) are coupled to r*(x)(crTcr^) (or T(x){o-^cr^i^}). 



Let us proceed with the analysis of the structure of the matrix 
V. Each of the four components of the matrix V, Eq. (l34l l. has 
the following block structure: 







(35) 



where the dimensions of the blocks (Vjkh, (V//t)L, and (Vjk)x 
are 6 X 9, 9 X 6, and 9x9, respectively. From Eqs. ([TJ, ( l34l l 
and ( l35T l, it is straightforward to obtain explicit expressions for 
the matrices (Vjk)a and (Vp^ (j k = 1,2, a l, x). We 
will instead introduce them implicitly - through their action on 
probe vectors. 

Let ai and 02 be 3-component column vectors. Using these 
vectors we create two 'probe' vectors 



fli <8i a2. 



fl2 



(36) 



The interaction matrices are now characterized by the following 
identities: 



(Vi2)-i(fll (802) 

(V2i)-i(ai (8)02) 

(yi2)-i(fli (» 02) 

(y2i)-i(fli 'Si 02) 

(Vi2)x(fll ® 02) 
(V2l)x(fll ® fl2) 

(yi2)x(ai ® 02) 
(y2*i)x(ai ® 02) 



(VI2). 
(V21), 

(K2)l 
(^2i)l 



fl2 

fll 

fl2 
fll 

fl2 
fll 

fl2 
fll 



2iTA<+>fl2(fli)2 

d 
d 

2iTA<+>a 1(02)2 

6 

-2/r*A<->ai(fl2)i 
-2iT*A<->a2(fli)i 

6 

-2rA<"'fli ® A(+>fl2, 
-2rA<+'fli ® A(">fl2, 
-2r*A<">fli ® A<+'fl2 
-2r'A(+>fli ® A<"'fl2 

ni «) (2;TA<+'fl2), 
2iTA<+'fli (g)ni, 
(-2iT*A*"^fli)®n2, 
in ® (-2iT*A<"'fl2), 



(37a) 
(37b) 
(37c) 

(37d) 

(37e) 
(37f) 
(37g) 
(37h) 

(37i) 
(37j) 
(37k) 
(371) 



where (oi), and (02), (and similarly for other three-component 
vectors in the following) refer to the i-th component (/ - 1,2,3) 
in the basis given by the choice of the single-atom Bloch vector 
(((T"), (cr+), (cr-')). Furthermore, the argument x in the coupling 
constants is for brevity dropped. Hi = (5,0, 0)^, «2 = (0, 5,0)^, 
d = (0,0,0)^, and 
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Note that the matrices V^, Vx, and transform the probe vec- 
tors, Eq. ( l36l l. to new vectors which can again be represented in 
the form of Eq. ( l36b . 



As regards the matrices A^"' and A*^\ they describe the cou- 
pHng of the Bloch vector of an individual atom to another atom 
- via coupling to the two-atom correlation functions, and the 
coupling of the latter back to the Bloch vector of the other atom. 
Remarkably, exactly the same matrices describe the coupling 
of a single atom to the classical laser field, see the components 
proportional to and O* in Eq. (|29] |. This is not a mere coinci- 
dence. We will prove below that, in the far-field limit, where the 
atoms exchange only single photon, the radiation of one atom 
onto the other atom is equivalent to that of a classical field with 
corresponding frequency distribution. 

3. The elastic component of CBS 

As obvious from Eqs. ( fT6b and (fTTT l. the elastic component 
of the spectrum results from the steady state of the two-atom 
Bloch vector, see Eq. (|6]i: 



{Q) = -(A + y)-'L+,,), 



(39) 



expanded in different powers of the dipole-dipole interaction V 
and finally configuration-averaged. 

3.1. Recurrence relations 

First, let us perform the expansion in powers of V. For this 
purpose, we split the 15 -dimensional vector (^) into a 6- and 
a 9- dimensional vector characterizing, respectively, the single- 
atom Bloch vectors and the correlations between both atoms: 



(40) 



where n corresponds to the power in the series expansion in V. 
By using the block structure of the matrix A + V and applying 
to it the formula for the inversion of block matrices |29] (see 
[Appendix A\ , one can show that these vectors obey the follow- 
ing system of recurrence relations: 



y 



(41a) 
(41b) 



with G+ = -M-\ Gx = -M-\ Va = (Vn)a + (V2i)« + (V^- + 
(y2j)a, a L, X. The initial condition reads: 

where Gj = The relations (HTI) represent a specific case 

of the general recurrence relations for the sub-blocks of the ma- 
trix (z - A - y)"', to be considered below in Sec. 14. II 

3.2. Steady state solutions 

Using the above recurrence relations of Sec. 13. II we will now 
explicitly calculate the steady state for fixed atomic coordinates, 
before we perform the configuration average in Sec. 13.31 

From Eqs. ( l40l i and (l42l) . we obtain the Bloch vectors for two 
non-interacting atoms (n = 0): 



The corresponding result for the correlations involves, accord- 
ing to Eq. (I41bl l. the matrix Gx, see also Eq. ( IA.3bl ). For the 
sake of reducing all calculations to the subspaces of atoms 1 
and 2, we shall use the following integral representation of the 



matrix Gx, proven in [Appendix B 
G 



n dcj' 
J-co 27r 



Gi(+/w')®G2(+'w'), 



(44) 



where Gi,2(±'w') = (+ia>' - Mi 2y \ The + (+) signs of the in- 
tegrands reflect invariance of the result, provided that the signs 
of the constituents are opposite. 

Application of the integral representation Eq. ( l44l i to 
Eq. ( IA.3bl ) yields: 

—Gi{i(jj')L®G2{-ioj')G2L 

CO 27T 

—Gx{ia)')Git®G2{-ioj')ll, (45) 
CO 27r 

Eq. (l45T l can be evaluated by means of the general sum rule 



valid for Re [z] - (see proof in Appendix C i: 



= ...Gi...G2(z).. 



(46) 



(0) 



(43) 



where . . . stand for arbitrary expressions which do not depend 
on (jj' and are identical in all three lines of Eq. (|46] |. Applying 
this sum rule to Eq. ( |45] | and using Eq. ( |43] | yields: 

{d-i®d-2f^ ^{d-if'^®{d-2f\ (47) 

which is the obvious result for two non-interacting atoms. 

The first-order correction to the vector ((T2)'^'* can be found 
from Eq. (HTt . Thereby, we arrive at 

<(?2)''* = 2/ri2(GiZ:)2G2A(+>G2Z^ 

-2iT2i(Gii)iG2A<"'G2i, (48) 

where Eqs. (I37ab . ( I37db . and W\\ have been used. The first- 
order correction for the Bloch vector of the other atom results 
accordingly as: 

<(?i)(i> =2jT2i(G2i;)2GiAWGiZ: 

-2;T*2(G2L)iGiA<-'GiL, (49) 

In order to be able to associate each term in Eqs. (l48l l and ( |49] l 
with the photon exchange processes depicted in Fig.lT] we have 
supplemented the coupling strength T with the corresponding 
indices, although their values are identical, i.e. T\2 - T2\ - T. 
As evident from Fig. [T] and in agreement with Eqs. (l48T l and 
$49[ . atom 2 is aff'ected by processes (a) and (d), correspond- 
ing to Ti2 and Tjj, whereas atom 1 is affected by (b) and (c), 
corresponding to T21 and Tjj, 

Similarly, we have derived the first- and second-order correc- 
tions to the correlation functions (<?! ® (T2)''^ and ((T2)*^^ (see 
Appendix D 1. Concerning the function (cri)*^', it follows from 
the expression (ID.3I I for the function (1x2)*^' after interchange 
of the lower indices 1 <-> 2 numbering the atoms. 



6 



3.3. Configuration average: ladder and crossed intensity 

Now, we have all ingredients to derive the elastic back- 
ground and interference contributions Lei and Cei to the average 
backscattered intensity given by Eqs. (fTSI l and ( fTTI l. Indeed, 
these ingredients, ((?i,2)*''\ (<?i,2)*'\ and (1x1,2)'^*, are given by 
Eqs. (|43]l, (|48]l, (l49]l and ( ID.3l l, respectively. We should finally 
perform configuration averaging of these expressions over the 
coordinates of both atoms. 

First, let us remind that both Lei and Cei are of second order 
in the interatomic couplings T{x) and T*{x). These couplings 
rapidly oscillate as a function of the distance - x/k between 
the atoms, see Eq. (|2]l. Hence, as already mentioned in Sec. 12.21 
only terms proportional to |r(ji:)|^ survive averaging over the 
interatomic distance, whereas all terms proportional to T(x)^ or 
T*(x)^ vanish. 

The next step is the average over the interatomic angle vari- 
ables. Also this averaging can be performed analytically, thanks 
to the following observation: the expression to be configura- 
tion averaged contains sums of products of matrix elements re- 



lated to either of the atoms. It turns out (see Appendix E 1 that 
all these matrix elements coincide, up to a position-dependent 
phase, with the ones obtained for an atom placed at the coordi- 
nate origin. This property reflects the fact that changing atomic 
positions changes nothing but the phases of the Rabi frequen- 
cies at the positions of atoms 1 and 2. The disorder average 
is then survived only by those terms who have not gotten the 
phase factor; all the rest terms vanish. 



Using the phase relations of Appendix E we can prove the 
following simple recipe in order to identify the terms surviv- 
ing the configuration average: we remind that, according to 
Eq. ( |34] |. T and T* indicate the exchange of negative and posi- 
tive frequency photons between the atoms. Moreover, the inten- 
sity of backscattered light, see Eq. (jSj, results from emission of 
a positive and a negative frequency photon from atom 1 or 2, in- 
dicated by the terms crj and o-^, j - 1 or 2, respectively. As we 
have found, exactly those processes survive the ensemble aver- 
age, where an atom emitting a positive (or negative) frequency 
photon has received this positive (or negative) frequency pho- 
ton from the other atom. For this reason, we have labelled in 
Sec. I3.2l the coupling strengths T12, T^^, T21 and Ljj in corre- 
spondence with the four processes depicted in Fig.[T] In particu- 
lar, the ladder component, which results, according to Eq. ( fTSI l 
from emission of a positive and a negative frequency photon 
from atom 2, involves only processes (a) and (d) of Fig. [T] The 
disorder average is hence performed by extracting, among all 
terms obtained when inserting Eqs. ( l43l l. (l48l l and (ID. 3b into 
Eq. (fTSI l. those which are proportional to ri2r2j. The crossed 
contribution, involving emission of a positive frequency pho- 
ton from atom 1 and a negative frequency photon from atom 2, 
consists of all the terms proportional to T^T^^, corresponding 
to processes (a) and (b) in Fig.[T] 

Introducing notations: 

(A(r(+;w))<^' = G{±iaj)h!~^^GL (50a) 
(A(r(w))<2' = GA<+'G(-/w)A<"'GL 

+ GA*"'G(/w)A<+>GL, (50b) 



and dropping the common prefactor 8ri2L|j, we arrive at the 
following final expression for the elastic background intensity 

Le, = {a-f\af\{Acr\Q)f\ar 

H-<Acr-(0))<2'(cr+)(°> + (Acr+(0))<+'<Acr-(0)>^"' 
H-<Acr+(0))<"'(Acr-(0))<+') 
dcij' 



+ 



2n 



r"'(w')«Ao-+(w'))"<o-"> 



\(0) 



where 
with 



-H<Ao--(w'))*^'(o-+)<°)), 
P'^°\u) = {G{-iu)6f)2 + {G{ioj)df)u 



(51) 
(52) 



^(0) 



= -/A^+^GL -H «2 - (GL)iGL, (53a) 

= +iA^-'GL + ni- (GL)2GL, (53b) 

gives the inelastic spectrum of single-atom resonance fluo- 
rescence 1 30] known also as the MoUow triplet (see also 



Appendix F I. 

For the interference contribution, we obtain similarly: 

Cei = {cr^f\cr-f\Acr-{0)r\Acr^{0)r^ 
+ <cr+)(°><cr+)(0)<A(r-(0)>(->(A(r-(0))W 

27T 



■f 



{Acr+(-icj')y-\GA^^>G(iaj')df'')2 



+ <crn(°) r ^(GA^--'G(-ic^')df\{Acr-(iaj')r\ 
J-00 2n 



(54) 



4. Inelastic spectrum of CBS 



4.1. General remarks 

In the previous section, we have presented a detailed deriva- 
tion of the elastic spectrum of CBS backscattered from two 
laser-driven atoms. The analytical results, Eqs. ( fSTT ) and ( l54l l. 
for the elastic ladder and crossed component are manifestly 
single-atom expressions. While the physical setting leading to 
such expressions will be discussed below, here we will proceed 
with the inelastic spectrum. Although we wifl need essentiaUy 
the same tools as we used for the derivation of the elastic spec- 
trum, the pathway to the final results is much lengthier for the 
inelastic case. Henceforth, we will restrict ourselves to a sketch 
of the derivation. 

In Sec. [3] it was shown that the background and interference 
elastic intensities as expressed by Eqs. ( 116117b correspond to 
processes (a),(d) and (a),(b) of Fig. [T] respectively. The expres- 
sions given by Eqs. ( l25b and ( |26] | generalize Eqs. (fTST i and ( fTTb 
to account for inelastic scattering, whereas the excitation ex- 
change processes remain the same. In other words, the ladder 
and crossed inelastic spectra, Eqs. ( l25T l and ( l26b . are also bilin- 
ear forms in the matrices V12, V^j and V\2, Vjj, respectively. 

Hence, to obtain the inelastic spectrum of CBS, we will solve 
Eq. ( fT9] ) for j = 2 using Laplace transformation, with the initial 
conditions being given by Eq. ( l20l l. perturbatively to second 
order in the matrix V. 



1 



4.1.1. Perturbative expansion of the resolvent 
The perturbative solution of Eq. ( fT9b reads: 

Ajf (z) = /?(">(z)Ai-*f' (0) + R^^\zW^^\0) 

+ 7;(2)(2)Aif' >(0), (55) 

where R^'\z), I - 0, 1,2, represent subsequent terms in the se- 
ries expansion of the resolvent matrix 



R(z) = (z - A - y)-' 



(56) 



in V. Obviously, the only difference of the perturbative expan- 
sion of R{z) from that of the matrix {-A - Vy^ is that, in the 
latter case, the Green matrices are evaluated at z = 0, while in 
the former case at z = -iv. Accordingly, the four sub-blocks of 
the matrix Riz)'. 



-Lx - V, 



z - Mx - Vx 



satisfy the recurrence relations generalizing Eqs. (HTt : 

Rf^(z) = G,(z)V.Rf,-'\z), 



r\%) = G,(z)V,R^^-'\z), 
Rf^(z) = G^L^Rf^(z) + G^V^R 



+ G^V,r\"-'\z), 



1,2, 



with the initial condition 



<'(z) = G+(z), Rf^(z) = Gx(z), 



(57) 

(58a) 
(58b) 

(58c) 
(59) 



mdR 



(«) 

kl 



for n < 0. 



4.1.2. Initial conditions 

We will next show that the vectors of the initial conditions 
As^^""'(0) on which the sub-blocks of /^'■"•'(z) act have the same 
structure as the probe vectors introduced in Eq. ( |36] |. 

Indeed, from the definition ( l20t of A.?2(0), one obtains 



what leads to the following expressions 



A^f" (0) = {cr^ay 



(0) 



A^-f (0) = (cr^f^ - {cr^f\^} 



y2) 

\(0)//)\(2) 



,(1) 



(60) 



(61a) 



(61b) 



(61c) 



The vectors (^)'^"* are manifestly of the form (l36l l by virtue of 
the relations (l37l l. (1411 1. and (|46] |. The same holds for the vector 
(o"2 2)'"* which is defined as 



in) 



iA^-\<T2Y"^ 
((Tl (g) /A<-'(?2)<"' j 



(62) 



Hence, the relations (l37T i. (HTt . and ( fSST l, supplemented by the 
integral representation (l44l i together with the identities derived 
in Appendix C are sufficient to evaluate all terms in Eq. i55[ . 



4.1.3. Configuration averaging 

The configuration averaging is reduced to the same proce- 
dure as for the elastic intensities. By choosing terms in the 
ladder and crossed spectra that are bilinear forms in the ma- 
trices Vi2, Vji ^^'^ ^12' ^12' respectively, and keeping only 
terms proportional to \T\^, we automatically select the compo- 
nents surviving the disorder averaging. The averaged expres- 
sions for ([Asf \-/v)]i>eonf. and ([Asf (-'V)]4e'^"''''^)conf. then 
follow from [As'2^'(-;v)]i and [As2^\-iv)]4, respectively, after 
dropping the atomic indices. 

4.2. Inelastic ladder spectrum 

After configuration averaging of both sides of Eq. ( |55] |. we 
obtain 



<[Arf\-/y)]i)eo„f. = L^^-^\-iv) + L^'''\-iy) 
+ L^"-\-iy), 



where 



L^^-"-"\-iy) = <[7?P-«)(-/y)A,?f >(0)]l)eonf.. 

In accordance with Eq. (l25T l. the ladder spectrum reads 
2 Re 



(63) 



(64) 



i-inel(v) 



0^'°>i-iv) + L''''>i-iv) + L'"''>(-iv)). (65) 



r(i;i) 



The explicit expressions for the functions L*^'"'(-!v), 
L*''''(-iv), and l!-^'^\—iv) are given in Appendix G 



Before we write down the final form of the inelastic ladder 
spectrum, let us introduce a new spectral function 

p(2\w', v) = 2 Re [{G{-iv)[6f{oj') 

+ A^-^G{-iv + /w')A*+'G(-!v)df 

+ A^->G(-iv + iu)')6^^\iu)') 

+ A(+'G(-/v - ioj')A^-^G{-iv)6f 

+ AWG(-/v - ia)')6^^\-ia)')])X 



(66) 



By using the function P^^Ho)' , v) and the inelastic single-atom 
power spectrum P^^\a)'), see Eq. ( |52] |. we can write down the 
expression for the inelastic ladder spectrum in a compact form: 



iine.(v)=f r^p(0)(c.')P<^' 

2n In 



(dL»',V) 



+ ±(^+)(0)(^-)(0)p(2)(0,y) 

In 

+ — (Acr+(/v))*"*(Acr-(-/v))(+'P*">(v) 
In 

+ — <Acr+(/v))*+>(Acr-(-/v))<"'P*"*(v). 

2n 



(67) 



Not only has the expression for the inelastic background spec- 
trum, Eq. ( l67l l. become short thanks to the function P^^\cj' , v). 
The latter function has a clear physical interpretation as the 
cross section for a weak probe beam of frequency u' to be scat- 
tered into the frequency v by a strongly driven atom. While 
this issue will be discussed in more detail when introducing the 
pump-probe approach in Sec. |5]below, we will next derive the 
inelastic interference spectrum. 



4.3. Inelastic crossed spectrum 

Repeating the same steps as taken previously in the deriva- 
tion of the inelastic background spectrum, we obtain 

<[A*f >(-/v)]4e*-'''^)conf. = CP^«'(-/y) + C^'^'^C-iv) 

+ C<"'^^(-iV), (68) 

where 

C(2-«;«)(_,v) = <[7;P-«)(-/y)A,?f '(0)]4e*-''-)eo„f., (69) 
from which the inelastic spectrum reads 

Ci„el(v) = ^(C<2-<'>(-/y) + d''"(-/v) + C<"'2'(-/y)), (70) 

n 

where the explicit expressions for the functions C^^'"*(-/v), 
C^^''*(-'V), and C*"'^*(-/v) are presented in |Appendix G| 

The crossed spectrum can conveniently be represented with 
the help of the frequency correlation functions P*^*(w', v) and 
P^-\a)',v) = (pW(v,w'))*: 

p(+'(w',v) = (G(/v)A*+>G(/v - /w')(?r)2 
+ {G{iv)&^\iu>'))2 



+ {G(-iv + ioj')A''^>G(-iv)df\ 
+ (G(-iv + ico')d'-^\iio'))u (71) 



P"(w',v) = iGi-iv)A^-'G(-iv + ico')df\ 

)i 



+ {G(-iv)d[-\-iaj'))i 
+ {G(iv - ia)')A^-'>G{iv)6f\ 
+ (G(iv - ia)')6Y\-ioj'))2, (72) 

The physical significance of the functions P^^Ho)' ,v) and 
P^^^o)' ,v) as spectral correlation functions of a two-level sys- 
tem subjected to a bichromatic driving will be clarified below 
within the same setting as the function P^^\to' , v). To this end, 
they allow us to write down a concise expression for the crossed 
inelastic spectrum as: 



C'inel ( 



27r J_„ 



dto' 



P*+V,v)P"(v- w',v) 



+ — (cr+)("\Acr-(-/v)/"^P^+^(0, V) 
+ — <cr-)W(Acr+(iV))<+^P^"^(0, v). 

2n 



(73) 



5. Pump-probe approach 

In Sections [3] and |4] we derived the configuration-averaged 
double scattering elastic and inelastic spectra of the background 
and interference contributions to CBS from two atoms driven 
by a laser field of arbitrary intensity. These quantities, given 
by Eqs. (ISTT i. ( l54b . d&Tl i. and ( l73T l. are expressed in terms 
of single-atom averages and spectral correlation functions of 
which some, e.g. pW(w',v), P^-\(jj',v), and P^^\u' ,v) have 
been introduced rather formally. We will next furnish them with 



a physical meaning. This meaning can be straightforwardly 
clarified because precisely the same expressions have been ob- 
tained by solving optical Bloch-lik e eq uations for a single atom 
subject to bichromatic driving ll2lll22ll . 

Here, the driving field £(f) is composed of the injected laser 
field and the photons scattered by the other atom: 



£(f) = Bie'"^'-' + V 



+ &le' 



+ 'e 



(74) 



both represented as a classical field with positive and negative 
frequency components. The laser field amplitude &i is pro- 
portional to the Rabi frequency Q. In correspondence to the 
approximation of large distance between the atoms, the probe 
fields v*^^ and v*"' are assumed to be much weaker than the laser 
field &L- As akeady argued in the introduction, this approxima- 
tion enables the modeling of the atomic radiation as a classical 
field acting on the other atom. We then consider the frequency 
correlation function 



C{0J\,(jD2) 



dtidt2 
In " 



X i{cT^(ti)cT-(t2)) - {crUti)}{(T-(t2))) , (75) 

representing a positive frequency photon oji and a negative fre- 
quency photon (jj2 (both with respect to the laser frequency) 
emitted by a single atom subject to the bichromatic driving, 
Eq. ( l74b . and expand it in powers of the weak probe field am- 
plitudes v^^' and v^"^: 



C|„..,=o = 5(wi-W2)^(°'(wi), 



dC 



dC 



- 6{(jJ2 - (jJ\ - tij)P^^\aj, 0)2), 

- 6{(X>\ - (X>2- LL>)P^ \tL>, tL>\), 



= 5(W1 - U>2)P^'^\u),U)i) 



(76a) 
(76b) 

(76c) 

(76d) 



,.(±1=0 



The 5- functions in Eqs. ( I76all76dl) originate from integrating 
over f+ :- t\ + t2 in Eq. (fTSl ). and are thus a consequence of 
time translation invariance or energy conservation. As we have 
found, the quantities defined by Eqs. ( l76aH76dl i - i.e. by the 
perturbative solution of the single-atom Boch equations with 
bichromatic driving - exactly correspond to the quantities f '"^ 
f (=t^ and f"^) introduced above in Eqs. ( I52I66I71I72| |. This es- 
tablishes the equivalence between the master equation and the 
pump-probe approach. 

Understanding the physical meaning of the quantities P^^^, 
f and f becomes easier with the help of diagrammatic lan- 
guage. Let us address a modified version of Fig. [1] plotted in 
Fig. 121 Namely, we added in Fig. |2] the laser field with which 
both atoms continuously interact, the positive and negative fre- 
quency components of the backscattered field, and combined 
the exchange of excitation diagrams (a),(d) and (a), (b) of Fig.[T] 
into the processes (a) and (b) of Fig. |2l respectively. As dis- 
cussed in Sec. 13.31 these are the processes which survive the 
configuration average, after solving the two-atom master equa- 
tion ([1) perturbatively to second order in the interaction matrix 
V. 
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t "I 



(a) 



(b) 



Figure 2: Double scattering processes contributing to (a) ladder and (b) crossed 
spectra of CBS. Thick an'ows indicate (positive and negative components of) 
the laser field driving the atoms (black spots), thin solid (dashed) arrows in- 
dicate positive (negative) frequency components of the fields scattered by the 
atoms. In all expressions, the frequencies v, ai' and ai" are relative with respect 
to the laser frequency w/,. 




(a-1) 



(a-2) 



the corrections (A(t(+/w))*''* and (A(?(<y))^^^ to the Bloch vec- 
tor defined in Eqs. dSOalSObb . This yields the other terms of the 
inelastic ladder component, Eq. (|67] |. where one atom scatters 
inelastically and the other one elastically. 

In the same spirit, the crossed spectrum is obtained as fol- 
lows: diagram (b-1) describes a positive frequency component 
of the field backscattered by the atom subjected to the laser 
field and to a positive frequency component at the frequency 
u)" emitted by the other atom. The atom also emits a nega- 
tive frequency component at the frequency w' towards another 
atom. The frequency distribution of these photons is described 
by the function P*^*(v, w')- The frequency w" is then deter- 
mined, due to energy conservation in the scattering processes, 
as (jj" - V - (1)' . Hence, diagram (b-2) is described by the 
function P'^ Hy - (J ,v). Integrating over w' yields the first 
term of Eq. (|73] ). corresponding to inelastic-inelastic scattering, 
whereas the remaining terms are obtained by including a single 
elastic event. 



6. Conclusion 



w 



(b-1) 



(b-2) 



Figure 3: Building blocks (diagrams) including the individual atoms with their 
incoming and outgoing field. The combined diagrams (a-1) and (a-2) give the 
ladder contribution (diagram (a) in Fig.|2), while the combined (b-1) and (b-2) 
diagrams give the crossed contribution (diagram (b) in Fig.|2). 



Now, let us split each of the diagrams (a) and (b) in Fig. [2] 
into 'building blocks' , in order to represent the individual atoms 
together with their incoming and outgoing fields. The resulting 
4 diagrams are shown in Fig. [3] 

Diagram (a-1) shows an atom subjected to a laser field. Con- 
sequently, its outgoing fields correspond to resonance fluores- 
cence with the frequency distribution given by the function 
P^^\ci)') defined in Eq. (l52b . On diagram (a-2), an atom is sub- 
jected to two classical fields: the laser field and a weak field at 
frequency (J scattered by the other atom. The response of the 
atom to such a bichromatic field expanded to second order in 
the weak field amplitude is described by the function f ^^'(w', v) 
introduced in Eq. ( |66] |. Integrating over the frequency w' leads 
to j'c/w'P<°*(w')^'*^*(<y', v), which is the first contribution to the 
inelastic ladder intensity, Eq. ( l67b . This is the contribution 
which originates when both atoms scatter inelastically, since the 
single-atom response as defined in Eq. ( fTsT l above includes only 
inelastic scattering. However, the same analysis can be repeated 
including also elastic scattering, i.e. considering Eq. ( fTSl ) with- 
out subtracting the term (cr^^')(cr^"') on the right hand side. As 
it turns out, the elastic response of a single atom to a probe field 
of frequency oj is then determined, in first and second order, by 



We have studied the spectrum of laser light doubly scattered 
between two scalar atoms. By solving the two-atom master 
equation to second order in the inverse distance between both 
atoms and averaging over the random positions of the atoms, we 
derived analytical expressions for the ladder (background) and 
crossed (interference) components of the coherent backscatter- 
ing signal. The results thereby obtained reduce to single-atom 
expressions coinciding exactly with the results based on the 
pump-probe approach to coherent backscattering |21, 221 ■ In 
the latter, heuristic, approach, all correlation functions are cal- 
culated by considering from the very beginning a single atom 
subject to a bichromatic classical driving field. 

This does not contradict the well-known fact (see, for in- 
stance, is, 28]) that atomic radiation exhibits photon antibunch- 
ing and, hence, is nonclassical. Photon antibunching and other 
nonclassical properties of resonance fluorescence are quantum 
statistical properties which do not reveal themselves inasmuch 
as single excitation exchange is involved. Instead, these prop- 
erties would become relevant for smaller interatomic distances, 
where several photons are exchanged between both atoms. In 
the case of two distant atoms, corresponding to the experimen- 
tally realistic case of a dilute atomic medium, however, the 
classical modeling is perfectly adequate - at least in order to 
describe the intensity of the atomic radiation. Whether an ap- 
propriately extended pump-probe approach could also be able 
to reproduce intensity-intensity correlations 1 2^ 31 1 is a subject 
for future studies. 

As shown in this paper, the classical modeling of the fields 
exchanged between the atoms allows to reduce quantum me- 
chanical calculations from the two-atom Hilbert space to single- 
atom Hilbert spaces. In combination with multiple scatter- 
ing theory as already existing for nonlinear classical scatterers 
[ 21, 151, the pump-probe approach thus opens the way to treat 
coherent backscattering of intense laser Ught by a dilute cloud 
consisting of a very large number of atoms. 
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Appendix A. Inversion of A + V 

Given a square matrix B, one can find its inverse by partition- 
ing it into 4 blocks, the diagonal blocks being square matrices, 
and applying the formula [Ref. L29.1 . par. 0.7.3] 



B 



[b,-b2b-^'bi\' b-'b2[b3b-'b2-b4] 
[b3b-'b2 - b4l' b^b-' [^4 - bib-'b2\ 



(A.l) 

-I \ 



For B = A + y. 



bi = M+, 
b2 = Vn, 

bi^L^ + y„ 

^4 = Mx + Vx- 



(A.2a) 
(A.2b) 
(A.2c) 
(A.2d) 



Expanding blocks B^\ and B^^ to second order in V, we estab- 
Ush the following relations 



(0) 



G2L 
Git 



<(?I®<?2)(0)=GxL: 



(<?2> 



(0) 



<C?2> 



(<?l) , 

(0) 



(<?2) 
(<?l) 



(A.3a) 
(A.3b) 
(A.3c) 

(0) 



(I) 



<C?2> 



(2) 



+ GxL^ 



(<?2) 
(<?l) 

<'?2) 



(2) 



(A.3d) 
(A.3e) 

+ Gxyx<(?i ® (?2>'" 

(A.3f) 



where G+ = -M;' = G2 ® Gi, Gx = -Mx', with G; = -MjK 
The recurrence relations ( IA.3b can be generalized to the form 
(HTt by the method of induction. 

Appendix B. Integral representation of the Green's matrix 
Gx(z) 



Since Mx = Mi ® / + / (g) M2, we can write 

„Mx/ _ „M,t ^ Mil 



Laplace transformation applied to the left and right hand side 
of Eq. (IB. lb gives 



J _ r 

:-Mx~ Jo 



dte-^'e^'' ®e^-'. Re [z] > (B.2) 



In order to transform the right-hand side of Eq. ( IB.2b to tensor 
products of Laplace transforms of individual atoms, we insert 
5{t - f') = 2b II doj'e'"'^'-''^ = ^ £^ doj'e''^'^''-'^ to get: 

Gx(z)=^J fl'w'Gi(|+/w')®G2(|+/w'), (B.3) 
with 

<^x(^) = ^ir' _ . (B.4) 

z-Mx \2 / i + iw'-Mj 

For z = 0, we obtain 

^Gi(+/w')®G2(W), (B.5) 
00 27r 

which coincides with Eq. (l44b . 

When there are two Green's matrices, the generalization of 
Eq. (ES) will be: 

...Gx(z)...Gx(z)...= 

J-co J-00 

...G,(|±/«')...G,(|±/..")... 
® . . . G2 + iw' j . . . G2 + ;W j 

Appendix C. Sum rule (l46t and related identities 

The eigenvalues /I/; of the Bloch matrix Mj given by Eq. ( |29] ) 
are defined from the characteristic polynomial [30.1 



f{z) = (z + 2r) ((z + yf + 6'-) + (z + r)"' 



(C.l) 



We note that Ak are independent of the atomic position, and 
Re [Ak\ < 0. 

Solving the eigenvalue equation, one finds the right (|m^» and 
left ((v^l) eigenvectors of the matrix Mj (which coincide with 
the right eigenvectors of the matrix Mj). Since (v^|m^) - for 
k + q (see [29], par.L4.6), one can represent Mj as 



where 



M 



(C.2) 



(C.3) 



(B.l) 



is the projector on the subspace corresponding to the eigenvalue 
Ak- One can easily deduce: 

{.z-Mj)^Y}^-A,)Pi, Gj(z)^J]jZJ-, ^^-^^ 
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With these preHminaries, we write down Eq. ( |46] | using 
Eq. (IC4l) 



Appendix D. Expressions for the functions ® (?2>*^' and 



J_oo In ^ a>' - ij + lAk 



X . . 



p2 



y 

^ y ^1 

J-oo 27r 4^ (w' - if + 



X . . 



OJ' + if - i/lfl 



(C.5) 



These integrals can be easily calculated using the residues the- 
orem. The integrands in Eq. ( IC.5b as a function of the com- 
plex variable to' have two simple poles at cu[ - iz/2 - iA^ and 
= -iz/2 + iAq. Taking into account that Re [z] = and 
that Re [Au] < we note that Im[w'[] > and ^[wj] < 0. 
By closing the integration contour in the upper half plane, we 
arrive after simple transformations at 



p2 

V '? 

'Z^7-A 



■G2(Z). 



(C.6) 



which completes the proof of Eq. ( l46b . 

As a concomitant result that we have used in our derivations, 
we will present here the following identity: 

... 02^-1^ -ico"y..G2(-ico')... (C.7) 



n dio' 



. Gi(iw') . . . G2(-!W - idLt') . . . G2(-i<^') . . . 



The proof of Eq. (IC.7b can be easily obtained by using the spec- 
tral decomposition, Eq. ( IC.4l i. and direct evaluation of the inte- 
gral over Lo" with the help of the residues theorem. 



In Appendix G there appears a number of integrals which 
can be calculated by using a slight modification of the sum rule 



dco' 
'2^ 



.G(zi)G(icj')...G(-icj'). 



XCO J / 
-^...G(ico')...Giz2)Gi-icj')... 
CO 2n 

= ...G(zi)...G(z2)..., (C.8) 

where zi and Z2 are independent from co' complex numbers with 
Re[zi] = Re[z2] =0. 



Using Eqs. ( I41I47I48I49T |, we obtain 

<(?! (g) (T2)*" =-2iT2*i(GiL)iGiL®G2A<">G2L 
+ 2iTi2(Git)2G 1 1 ® G2 A(+>G2Z! 
+ 2iT2i (G2L)2G 1 A<+'G 1 L ® G2L 
-2iTi2(G2L)iGiA<"'GiL®G2L 

-2iT2^£^G,(iV)(?f, 

® G2(-iw')A*"*G2L 

+ 2iTn^ '^G,{-iJ)6% 

® G2(!w')A^^^G2L 

— Gi(iw')A' ^GiL 
00 2n 

® G2(-/w')G^°' 



2iT,*2 r 

%J —0(. 



)G2(-/«')G^";, 



2;2 

doj' 
'2^ 

m 



■Gi(-iio')A'-''GiL 



(D.l) 



where we have introduced new vectors 



C^fl = -iA*+>G/ + n2- {Gjt)xGjt, 
dfl = -H/A<->G;L + ni- (GjL)2GjL. 



(D.2a) 
(D.2b) 



The vectors G^?J and G^'j play an important role in our deriva- 
tions: they represent the initial conditions for correlation func- 
tions used in the calculations of the inelastic spectrum in single- 



atom resonance fluorescence (see Appendix F 1 



From Eqs. ( |46] |. ( ID. Il l, (ID. 2b . and (14 lab , the second-order cor- 
rection ((T2)'^* for atom 2 follows straightforwardly. Here, we 
write down detailed expressions only for the terms proportional 
to |rp, and do not specify those proportional to or (T*)^ 
which will vanish after performing the configuration average, 
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see Sec. 13.31 

{d-2f^ = + 4ri2r*i(GiL)2(GiL)iG2A(-)G2AWG2L (D.3) 
+ 4ri2r*i(GiL)2(GiL)iG2AWG2A"G2L 
+ ATnTUGit) 1 (G i A'-'G , z:)2G2 A'+'GzZ! 
+ 4r2ir*i(G2L)2(GiAWGiL)iG2A«G2L 

X G2A(">G2(/w')A^+>G2L 

X G2A<+)G2(-/w')A*"'G2L 

+ 471271*2^ ^(Gi(/d^')A(">GiL)2 

X G2A<+>G2(-/w')(?2;l 

^(G,(/d^')AWGiL)j 

X G2A<->G2(-/W')(?2;2 + (^^) ' " ' + ^T*f .... 

By reasons of symmetry, the corresponding result (cri )*^' for 
atom 1 is obtained by exchanging the Green's functions G\ «-> 
G2, the vectors Gj"j <-» Gj'j and Gj^j <-> Gj'^j' the coupling 
strengths T\2 <^ T21 and r*^ <-> Tjj. 

Appendix E. Phase relations 

The following relations between expressions for an atom 
placed at position r, (left-hand side) and an atom placed at the 
origin (right-hand side) can be obtained by direct calculation: 





(GjLh 


= (GL)ie*^''', 


(E.la) 




(GjL)2 


= (GL)2e"'*'''^ 


(E.lb) 




''Gjt), 


= (G(z)A(+>Gz!)ie2'*^-0, 


(E.lc) 




''GjL)2 


= (G(z)A(+>GL)2 


(E.ld) 


(G/z)A( 


-^GjL), 


= (G(z)A"GZ:)i, 


(E.le) 


(G/z)A( 


->G/)2 


= (G(z)A(->GL)2e"^*"'''^' 


(E.lf) 


(G;A(->G/z)A( 


^>G/)i 


= (GA(->G(z)AWGz!)i/'-''^ 


(E.lg) 


(G;A(->G/z)A( 


^>G/)2 


= (GA(->G(z)A<+'GL)2e"*"-''' 


(E.lh) 


(G,A(+>G/z)A( 


->g/)i 


= (GA(+>G(z)A<->GL)ie''*^'--''^ 


(E.li) 


(G,A(+>G/z)A( 


->G/)2 


= (GA*+*G(z)A("'GL)2e"'*"'''^ 


(E.lj) 



where G(z) = Gy(z)|^^g, and z is an imaginary number. 

Each of the terms to be averaged represents a product of 
matrix elements of the form given by (IE. Il l (see, for instance, 
Eq. dP.ll l). Then, using dE. II ). one can easily deduce the to- 
tal phase of each term. There are two types of terms surviving 
the disorder averaging: (i) terms without the total phase factor 
and (ii) terms with the total phase factor e^''^'-'''^- . The former 



contribute to the ladder spectrum, and the latter to the crossed 
spectrum. This rule can also be generalized for the inelastic 
spectrum. 



Appendix F. Inelastic spectrum in single-atom resonance 
fluorescence 



The expression for the inelastic spectrum, Eq. (l22l l. applied 
to the single-atom case can be written as 

1 r°° 

^<'"(v)=-J^ «'T(e'-(A.V(T))j (F.l) 
+e-"' {A?j.,-{T)\), (F.2) 

where Asj-^ij) obeys the equation of motion 

A^j.^ = MjAsj.^, (F.3) 

with initial conditions 

A^j-AO) = <cr}(?;) - (crtXc?,), (F.4a) 
A^j.^O) = (cfjcrj) - {cT-j){&j). (R4b) 

Comparing Eqs. (IF.4b and (ID. 2b one obtains 

A.?,;+(0) = (F.5) 
A,?,;_(0) = C?^°J. (F.6) 

Solving Eq. (IF.3I I with the help of the Laplace transform, we 
obtain the following expression for the inelastic spectrum of 
resonance fluorescence: 

P'^v) = ^jGjiiy)d%\ + (G,(-/v)d^°J)J. (R7) 

The spectral distribution of resonance fluorescence is indepen- 
dent of the coordinate of the atom. Accordingly, the atomic 
index j can be dropped in the right hand side of Eq. ( IF.7I ). since 

(g//v)G^)^ = (G(/y)(?f , (R8a) 
{Gj{-iv)df)^ = {G{-iv)6f\ , (R8b) 

where the expressions with the dropped atomic index on the 
right hand side of Eq. ( IF. 8b correspond to an atom placed at the 
coordinate origin. 



Appendix G. Spectral correlation functions L*^ "'"'(-iv) 
and C*2-n;n)(_jV) 



Omitting the common prefactor 4|rp, we get the following 
results for the functions L^^^"'-"\-iv): 



L<2;0)(-/v) = (GL)i(GL)2 (G(-/v)A<+>G(-/v)A<-^G(-/v)(?f )j 
+ {Gt\{Gt)2 (G(-/y)A(-^G2(-/v)A(+^G(-/v)(?f )j 

: (G(-/v)A(+>G(-/v - /w')A<"'G(-/v)(?f )j 

(G(-/v)A<->G(-/v+ /w')A(^'G(-/v)(?f )j , 



+ 
+ 

XI 



(G.l) 



+ 

XI 



L(i'i>(-/y) = (GL)i(GL)2 (G(-/v)A(+'G(-/y)(?^^(0))^ 
+ (GL)i(GL)2 (G(-/y)A<->G(-/v)(?^^'(0))j 

(G(-/v)A<+'G(-/v - /w')(?2"'(-iw'))j 

X (G(-iV)A(-'G(-/v + icjj')6\^\icjj'))^ 

:(G(-/w')A<"'GL)2(G(-/v)A<+>Gl)j 

X]t^(G(-W-/c.')GT), 
X (G(/w')A(+'GL)2(G(-/v)A(-'GL)i, (G.2) 



+ 

x| 

+ 



L(0;2)(-/v) = (GL)i(GL)2(G(-/v)C?f (0))j 

+ £ ^ (G(-/..')(?f (G(-/v)(?f (o.')), 
+ £ ^ (G(/c.')GT), (G(-/v)G^f (c'))^ 

X {G(-iv)G(-iLtj')A'-^Gtj (G(/w')A'+>Gz!) 

X (G(-;v)G(/w')A'^^GL) (G(-/w')A'"^GL) 

X (G(-/v)G(/w')A*+>GL) (G(-/w')A*">GL) 

X (G(-/v)G(-/w')A*"*Gl) (G(/w')A'^^Gl) , 

(G.3) 

where the vectors Gj '(-/v), G2^^(/v), and Gf\v) at w = are 

related to the vectors of initial conditions Ai^^'(O). They are 
defined as 

d^'\-iv) = /A<"'G(-;v)A(">GL-G(-/v)A<"'GL(gl)2 

-G£(G(-/v)A<-'GZ^)2, (G.4) 

d^^\iv) = /A("'G(/v)A(+>GL - (G(/v)A<+'GL)2 GL 

-[GL)^G(iv)A^^^GL, (G.5) 

^(v) = /A<"'GA*+*G(-/v)A<"'GL + /A<"'GA<">G(/v)A<+*GL 

- Gl(gA<+)G(-;V)A<-'Gl)2 - Gl(GA'-'G(/v)A*+'Gl)2 

- (GL)^ GA(+'G(-/v)A(->GL - (GL)^ GA<-'G(!v)A<+*GL 

- (G(/v)A<+'Gl)2 G(-/v)A<-^GL 

- (G(-/v)A<->GL)2 G(/v)A<+^GL. (G.6) 

To obtain the complex conjugated part of any of the expressions 
from Eqs. ( IG.llG.2IG.3IG.4IG.5IG.6l l needed for calculation of 
their real parts, one has to make the following interchanges in 
these equations: (a) / -/; (b) all subscripts 1 <-> 2; (c) all su- 
perscripts (+) <-> (-). The superscripts (0) and (2) should be left 
unchanged. These interchanges imply introducing 3 new vec- 
tors, G^^\iv), G^ \-iv), and G^^\v), whose explicit expres- 
sions follow immediately from Eq. (IG.4IG.5IG.6I ) by applying 
the substitutions (a)-(c) specified above. 

Next, we note that the 2 lower lines in Eq. ( IG.2I ) and 4 lower 
lines in Eq. ( IG.3I ) form, with their complex conjugates, pairs of 
terms for which integration over oj' can be performed analyti- 
cally with the help of Eqs. (l46llC.7IC.8T l. After these transfor- 
mations, the corresponding 12 integrals are reduced to 4 alge- 
braic expressions. 
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For the correlation functions entering Eq. (iTOl i we obtain 
C<2;0)(-/v) = (GL)2(G(-/v)A*-'Gl)^ (G(-/v)A*+)G(-/v)(?f )j 



+ £^(G(-/v)A(-)G(-/c.')GT)i 
X (G(-iV + !a;')A^^'G(-(v)df )^ , 



(G.7) 



C^'''\-iv) = (GL)2 (G(-/y)A(->GL)j (G(-/y)d^+'(0))^ 

+ £^(G(-/v)A(-)G(-/c.')GT), 
x(G(-/v + /w')(?2'"^(w')), 

^(GL),(G(-/y)A«G(-/v)G(-/«')(?f); 
X (G(iw')A*+'GL)2 , (G.8) 



C(0;2)( 



x| 
xl 



J-oo 27r 

x (G(-/v)G(-/w')A<"^Gl)j (G(/w')(?i''^(0))2 

X (G(-/v)G(-/w')(?2"^(0))j (G(/w')A<+>Gl)2 

X°° dco' r°° c/oj" 

(G(-/v)G(/w")A<->G(-/w')(?f )j 
(G(-/w")G(/w')A*+*Gl)2 (GL)i 

(G(-/v)G(/w")A*"'G(-/w')(?f )j 
c/w' r°° doj" 
oo 27r J_oo 27r 
(G(-/v)G(/w")(?^"'(-/w'))j 
(G(-/w")A<+>G(/w')(?r)2 

(GL)2(G(-/v)G(/w")G(-/w')A*"*Gl)j . (G.9) 



x| 
xl 



Calculation of the real part of Eqs. (IG.7IG.8IG.9b is done by 
employing the same substitutions as in the case of the back- 
ground spectrum (see above). The final result follows after 
subsequent integration of Eqs. (IG.7IG.8IG.9b over co" and, par- 
tially, over a>' by using the identities ( l46l l, ( IC.7b . and ( IC.Sb . 
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